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, We make theoretical predictions for the recently introduced variable (jf corresponding 

to the azimuthal angle between leptons produced in the Drell-Yan process at the 
LHC. As a consequence of this work we are also able to generate results for the more 
commonly studied transverse momentum Qt of the lepton pair. Comparisons of these 
purely perturbative estimates for the Qt case yield good agreement with ATLAS and 
CMS data, as we demonstrate. We anticipate that this work will help stimulate 
^ . measurements off/ at llie LUC. 
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Abstract 



1 Introduction 



The transverse momentum (Qt) distribution of lepton pairs produced in hadron-hadron collisions 
is one of the most studied observables in particle phenomenology. Perturbative QCD corrections 
are under control at next-to-leading order (NLO) and fully differential codes to generate them 
are publicly available, e.g. [THS]. However, fixed-order predictions fail to describe the small-Qr 
region and resummation techniques must be used in order to obtain reliable predictions. The 
resummation of large logarithms in the Qt spectrum has been studied by several groups [B1-H9j . 
The state of the art in the resummation is next-to-next-to-leading logarithmic accuracy (NNLL) 
in the resummed exponent |15| . 

On the experimental side, the transverse momentum distribution of Z bosons decaying to 
lepton pairs has been intensively studied by the and CDF collaborations at the Tevatron [20H22] . 
More recently, measurements have been performed at the LHC as well [231124] . 

State-of-the-art theoretical calculations have been compared to data coming from the Teva- 
tron experiments with the aim of testing QCD predictions and to assess the importance of non- 
perturbative physics. For instance, in Ref. [T7] a NNLL resummation, matched to NLO calculation 
was found to provide a good description of Tevatron Run II data, within theoretical and exper- 
imental uncertainties, with no apparent need for non-perturbative contributions. A somewhat 
different approach is used in the program RESBOS In this case the evaluation of the 

theoretical uncertainty that affects the resummation is less thorough and the agreement with the 
data is restored thanks to a non-perturbative contribution [22,26]. In our view, the difference 
between the pure perturbative result for the Qt spectrum and the one with non-perturbative 
corrections included can be similar in size to the perturbative uncertainty. We would thus prefer 
to evaluate the perturbative uncertainty carefully prior to including non-perturbative effects via 
a phenomenological study in conjunction with experimental data. 

Further obscuring the picture is the experimental resolution which affects the measurement of 
transverse momenta. For this reason novel variables have been introduced in [;27,28j and measured 
by the DO collaboration [29]. These variables, labelled ay and cj>*, both crucially depend on the 
azimuthal angle A(fr between the final state leptons, at low Qt- The experimental resolution for a*r 
and 4>* is significantly better than the one for Qt [28] . which enables the possibility to better test 
theoretical ideas and constrain non-perturbative effects. An example of such a phenomenological 
study with the RESBOS code can be found in Ref. [29J, where the data clearly disfavoured the 
inclusion of non-perturbative small- x broadening terms included in the RESBOS code. 

In previous papers |3UH32j we provided the details of a resummed treatment of these new 
variables and discussed their relationship to Qt and to each other. We computed the resummation 
to NNLL accuracy, carried out the matching to fixed-order NLO results from MCFM [1] and 
performed a complete phenomenological study, including a faithful estimate of the theoretical 
uncertainty. We compared our prediction to the data and we found that resummed perturbation 
theory provides an excellent description of the 4>* distribution with little variation in the quality 
of our description in the different rapidity regions. 

In this Letter we extend our work by providing resummed and matched predictions for the 
LHC. We begin by considering the standard Qt spectrum and we compare our results to the mea- 
surements performed by ATLAS and CMS. Having thus obtained a validation of our method, we 
then provide theoretical predictions for the </>* variable, which should stimulate its measurement 
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at the LHC. 



2 Resummation of the Qt and 0* distributions 

The resummation formalism for the Drell-Yan Qt spectrum is well established. The resummation 
can be performed to NNLL accuracy [15J and it can be matched to NLO calculations to obtain 
accurate predictions over a vast range of Qt- 

Here we merely remind the reader that in order to properly treat transverse momentum con- 
servation in the case of n-gluon emissions with transverse momentum k T ii one usually introduces a 
two-dimensional impact parameter 6, which is the Fourier conjugate of the lepton-pair transverse 
momentum Q T . The delta function which enforces momentum conservation in the transverse 
plane can be then expressed in a factorised form: 



i \ n 

Zh Ti + Q T )= A / d 2 be^r H^-kTi. (1) 
=i / 11 J i=i 



If one is interested in the Qt spectrum one integrates over the angle between b and Q T , obtaining 
a Bessel function. The resummed differential distribution has the form 

-777- (Qt, M, cos 9*, y) = ^ f° dbbQ T J Q (bQ T ) e -*(MW«) 
g\Qt sN c J 

x S (xi, x 2 , cos 6*,b, M, hq,hr, Hf) , (2) 



where x\ 2 = M=e ±y , being M and y the dilepton's invariant mass and rapidity, respectively, 

and b = Hp-. The resummed result depends on three arbitrary scales: renormalisation and 
factorisation scales, and fiF, as well as the resummation scale /xq. Variations of these scales 
around the dilepton invariant mass M provide us with an estimate of the theoretical uncertainty. 
In Eq. ([2]) the dependence upon the large logarithms we wish to resum is encoded in the radiator: 

flf&^A^oW) = L 9 (iHa s L) +9 w(* s LA,^)+%®(a s L,M,^), 

(3) 

where L = ln(b 2 [iQ) and a s = a s {nR), as well as in the cross-section E, via DGLAP evolution 
of the parton distribution functions from the hard scale of the process to Q = J~ . We refer the 

reader to the Appendices of our previous work [32] for the explicit expressions of the functions 
as well as for the cross-section X . It is well known that the 6-integral in Eq. ([5]) has issues 

both at large and small b. We evaluate the integral by introducing an upper cut-off 6 max and by 

freezing the radiator for b < 6 m i n , as explained in detail in Ref. [52] , 

An issue that has become apparent in the current study is related to the behaviour of the 

parton distribution functions at low momentum or, equivalently, at large b. In our Tevatron 

study [32J we decided to freeze the parton densities below Qq = 1 GeV. We tested this prescription 



1 Om expression for S does not include the O (a,) correction which is in principle needed to achieve complete 
NNLL accuracy. 
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by varying the value of Qo within a factor of two, finding that the sensitivity was much less than 
the perturbative uncertainty. However, at the LHC energies, the parton distribution functions 
fi(x,Q) are typically probed at lower values of x, where the Q dependence is steeper. In the 



for low values of the ratio fip/fj,Q the introduction of an abrupt freezing point at 1 GeV gives 
rise to an oscillatory behaviour of the Qr-distribution. We circumvent this issue by constructing 
an exponential extrapolation of the parton densities below 1 GeV. This procedure cures the 
oscillations but still produces curves with a peak which is noticeably shifted to the right with 
respect to our central value /j,q = fip = M-R = M. Such curves that critically depend on the 
behaviour of the parton densities below 1 GeV are not reliable in our framework based on collinear 
factorisation and DGLAP evolution of collinear parton distribution functions. It is far from 
obvious that one can trust this framework at very low scales, where small- x effects should be 
taken into account and collinear factorisation might not be valid, signalling the need for transverse 
momentum dependent parton distributions [33-35J. It would be interesting to see whether the 
sensitivity we observe, when we vary all scales independently, is also present in different theoretical 
implementations. However, only the approach of Refs |12H17| takes into account independent scale 
variation, as we do, and, to our knowledge, results at LHC energies have been published by that 
group for the Higgs Qt distribution but not for the Z transverse momentum. 

In this study we adopt a more phenomenological viewpoint. We note that if we vary the three 
scales (hq,hr and fxp) independently, all the curves for which ^ > 1 group together around the 

central one, while the ones for which — = i tend to form a distinctly outlying family of curves, 
with a shifted peak. Therefore, we decide to vary the perturbative scales independently, with 
the additional constraint ^ > 1, taking the resulting band as an estimate of the perturbative 
uncertainty of our resummed prediction. Hence the curves we produce are indeed insensitive to 
the behaviour of the parton distributions at very low scales. 

In this Letter we are also considering the variable eft*, which is defined as |28j 



where Acj) is the azimuthal angle between the two leptons produced by the Z/j* decay and sin#* 
is the scattering angle of the dileptons with respect to the beam, in the boosted frame where the 
leptons are aligned. This definition of 9* avoids the necessity of measuring magnitudes of lepton 
momenta and significantly helps the experimental resolution as explained in Ref. [28J. We have 
also introduced the angles <j>i between the momentum of gluon i with respect to the lepton axis 
in the transverse plane. Thus, <f>* is essentially determined by one component of the transverse 
vector Q T . For this reason, the formalism of Qt resummation can be applied to the variable <jf 
as well [30-32J. We can straightforwardly perform the integral over one component of b, obtaining 
a cosine instead of the Bessel function 



resummed calculation the argument of the parton densities is Q 



= we have noticed that 





x £ (xi,X2,cos9*,b,M,fiQ,fi R ,fj, F ) . 



(5) 
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The theoretical prediction for the Qt or <f>* distributions is obtained by matching the re- 
summed calculation to a fixed order computed with the program MCFM [I]: 

V°V matched V°WnNLL V°WnLO V<W expanded' 

with v = Qt,4>*- The last term in the above equation is the expansion of the resummation to 
O (a 2 ) and avoids double counting. We note that the difference between the resummation and its 
expansion is well behaved only up to values of Qt of the order of the Z mass. This can be linked 
to the fact that resummation formalisms in impact parameter space suffer from a non-physical 
behaviour at sufficiently large transverse momenta [17J. Therefore, for Qt > mz our result is the 
pure NLO. The expressions in Eqs. ([2]), ([5]), as well as the fixed-order prediction from MCFM, are 
fully differential in the dilepton kinematics, so that we can take into account any experimental 
cuts. 



3 The Q T and 0* distributions at the LHC 

In this section we compare our resummed and matched results to the measurement of the Z/j* 
transverse momentum distribution in proton-proton collisions at 7 TeV, performed by the ATLAS 
and CMS collaborations. The cuts applied by the two experiments slightly differ. We stress once 
more that because our predictions are fully differential in the lepton momenta we can take these 
cuts into account, without any need to extrapolate to the full lepton acceptance. Performing 
comparisons to the data in the fiducial region is very important if one wants to fully exploit the 
potential of the measurement, as recently pointed out, for instance, in Ref [I] in the context of 
production and studies of parton distribution functions. 

In particular, we compare our theoretical prediction to the ATLAS data (electrons and muons) 
in the fiducial volume, defined by p l T > 20 GeV and |r/| < 2.4, in the lepton-pair invariant mass 
window 66 < ran < 116 GeV. The CMS data instead are presented for p l T > 20 GeV, \r] l \ < 2.1 
and 60 < mu < 120 GeV. In Figure [1] we compare our theoretical prediction to the ATLAS and 
CMS data. The theoretical curve is computed using the CTEQ6m set of parton densities [36], with 
the value of the strong coupling taken from the fit, a s (Mz) = 0.1179. The curve is normalised to 
its own area. The band represents our assessment of the perturbative uncertainty. It is obtained 
by varying the perturbative scales /iQ, [Xf and between M/2 and 2M, with the constraints 
^ < — < 2 and — > 1, where i,j € {F,Q,R}. From this we can estimate the perturbative 
uncertainty to be around 10% in the peak region we are most interested in. Another source of 
theoretical uncertainty comes from the parton distribution functions. We estimate this one to be 
smaller (of the order 2% at low Qt) than the one obtained with scale variation. This is expected 
because we are considering the shape of the Qt distribution and the uncertainty coming from 
the parton densities largely cancels when we perform the ratio to the inclusive cross-section. Our 
predictions are computed without any explicit non-perturbative effects. We have checked that 
the inclusion of a Gaussian contribution in 6-space with coefficient c?np = 0.5 GeV 2 , taken from 
our Tevatron study [52] . produces distributions compatible with our perturbative uncertainty. 

The ratio plots show that we have a good description of the data in the Iow-Qt region, within 
theoretical and experimental uncertainty. For intermediate values of the transverse momentum 
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Figure 1: Comparison of the theoretical prediction NNLL+NLO for the Qt spectrum to the 
experimental data collected by the ATLAS collaboration (on the left) and by the CMS collabo- 
ration (on the right). The theoretical uncertainty is obtained by varying the perturbative scales, 
as explained in the text. The curves are normalised to unit the data are. 
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our prediction slightly overshoots the data points from both experiments, by 5% to 10%. In the 
large Qt region, where our prediction reduces to NLO, we correctly describe the ATLAS data. 
The agreement with the CMS data is instead not as good in this region; this feature is also present 
in the comparison to the NLO distribution performed by the collaboration itself [24J. 

The experimental collaborations compared their data [23, 24J to the fixed-order prediction 
computed with FEWZ [2] and to various Monte Carlo parton showers with or without inter- 
faces with NLO codes, e.g. POWHEG |37| . As expected, fixed-order perturbation theory works 
reasonably well at high-Qj- but resummation is needed in the Iow-Qt region. Parton showers 
and RESBOS [25] were found to provide very good descriptions of the data, although the main 
source of theoretical uncertainty at Iow-Qt, the resummation scale [Aq, was not explored in those 
studies. Furthermore, both parton showers and RESBOS predictions contain non-perturbative 
effects. 

Finally, a comparison between the theory prediction using SCET and the ATLAS data have 
been performed in |19j in the intermediate-to-low-Qy region. Good agreement was found in the 
study with the introduction of a non-perturbative contribution to describe the lowest-Qr data 
points. 

Having validated our resummation procedure at the LHC energies, we can now provide a 
prediction for the (/)* distribution in proton-proton collision at 7 TeV. We nominally choose the 
same selection cuts as the ones adopted by the ATLAS and CMS collaborations for Qt- The 
resummed and matched results are shown in Fig. [2J in the case of ATLAS cuts (on the left) and 
CMS cuts (on the right). The theoretical uncertainty is estimated as explained above, and in the 
low-0* region it is 0(10%). The size of the uncertainty band is comparable to the one we found for 
this distribution at the Tevatron. Furthermore, the height of the plateau is significantly reduced. 
This is due to fact that more radiation is produced when the parton distribution functions are 
probed at lower values of x, with the result that configurations with Born kinematics are less 
likely. 

We advocate the measurement of the <f>* distribution at the LHC, in order to accurately probe 
the \cw-Qt domain of lepton pairs produced via the Drell-Yan mechanism. Moreover, the LHC 
experiments can probe different regions of phase-space such as, for instance, forward rapidities 
or low values of the dilepton invariant mass. In these regions one expects small- x contributions, 
which go beyond the resummation formalism currently used, to become important, signalling the 
need for better theoretical modelling |38H41j . 

4 Conclusions 

In this Letter we have compared a resummed and matched prediction for the Qt spectrum of 
Z boson in proton-proton collisions at 7 TeV. The accuracy of the resummation is NNLL and 
it is matched to a NLO calculation obtained with the program MCFM pQ. Our prediction is 
fully differential in the leptons' momenta and therefore can be compared directly to the data 
in the fiducial region. We have found good agreement between the experimental data and our 
calculation in the 1ow-Qt region we are most interested in even without direct inclusion of non- 
perturbative effects, which is a similar finding to that of our previous phenomenological study of 
<f>* at the Tevatron [32] , Additionally in this Letter we have provided the theoretical prediction 
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Figure 2: Theoretical predictions at NNLL+NLO for the normalised (ft* distribution in pp colli- 
sions at 7 TeV. The kinematical cuts are the same as the ones adopted for the Qt case by the 
ATLAS collaboration (on the left) and by the CMS collaboration (on the right). 



for the (ft* variable at the LHC. 

Our resummed predictions are obtained with a computer program that we plan to release soon. 
The code can be easily interfaced with a fixed-order program (MCFM in this case), which provides 
the Born- level results. Our program then computes the desired observable in terms of the final 
states' kinematics and re-weights the Born cross-section with the resummation. For this reason 
the numerical code can be easily applied to different processes, which involve a colour singlet in 
the final state. For instance, the resummation for separation in azimuth of the final state leptons 
A(ft can be immediately obtained. Moreover, because the resummation is blind to the actual 
nature of the (colour-singlet) final state, the same code can be applied to compute resummed 
prediction for processes with different final states, e.g. the azimuthal separation between vector 
bosons or Higgs and Z bosons in associated Higgs production. The inclusion of gluon-induced 
processes, e.g. Higgs Qr-spectrum or azimuthal diphoton distributions, is also possible but it 
requires modifications of the resummed exponent. However, new structures arise if one wants to 
go beyond the leading-logarithmic accuracy [32] . Finally, the methods we have used here are also 
applicable to the study of azimuthal angle between coloured particles in the final state (jets in 
dijet production) which would also be interesting to explore in detail phenomenologically [43J. 
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